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QCD coupling up to third order in standard and analytic 

perturbation theories 



B.A. Magradze 



1 Introduction 



In publications |fl, |2| a new approach to perturbative QCD, the renormalization group (RG) invariant 
analytic approach (IAA) was proposed (for the exhaustive review on this approach we recommend 
paper ||). This method consistently takes into account the renormalization invariance and analyt- 
icity in perturbation theory (PT). 

In work a particular version of IAA has been formulated. According this version, analytic 
perturbation theory (APT), the observables in the space- like region are represented as a non-power 
, series in the special universal functions A n (Q 2 , f) (n=l,2. . . and f denotes number of quark flavor) 
||]. Analogical set of functions {2l n (s, /)}, s > 0, was introduced in the time-like region Q. Both 
sets of functions are determined by the QCD running coupling a s (Q 2 ,f) and can be calculated in 
APT analytically or numerically. A systematic mathematical investigation of these functions have 
been undertaken in works [5-7], in particular the oscillating behavior in the infrared region was 
■^j- ■ established. To calculate these functions beyond the one-loop level the iterative approximation for 
the coupling was used [||, ||, ||], or RG equation for the coupling has been solved numerically in the 
complex domain 



In papers [11, 12] the RG equation for the QCD running coupling, at the two-loop order, was 
solved explicitly as a function of the scale. The solution was written in terms of the Lambert W 
function. The three-loop order solution (with Pade transformed /5-function) also was expressed in 
terms of the same function [jl2] ]. Using the explicit two-loop solution, the analytical structure of 
the coupling in the complex Q 2 -plane has been determined [11, [j~2]| . The analytical formulae for the 



Oh, corresponding spectral function was found. Then the analytically improved coupling Jl^, [l4| was 
Q_i! reconstructed. 

Afterwards, in paper [15], the running coupling of an arbitrary renormalization scheme, to the 



k-th order (k > 3), was expanded as a power series in the scheme independent explicit two-loop 
order solution. The new method for reducing the scheme ambiguity for the QCD observables has 
been proposed in this work. A similar expansion for the single scale dependent observable, motivated 
differently, has been suggested in |l(| . 

In Sec. 2 and Sec. 3 we use the explicit solutions for the running coupling to calculate the universal 
quantities A n (Q 2 ,f) and 2t n (s,/) beyond the one loop order. The results for %l n (s,f), to second 
and third orders, are presented in the analytical form. 

In Sec. 4, the matching conditions for crossing the quark flavor thresholds for the MS scheme 
running coupling a s (Q, /), to the three loops, are solved analytically. By the way we construct the 
global (independent on f) universal functions, A n (Q 2 ) and 2l n (s) (both introduced in ||). These 
functions can be used in the whole momentum space. 

In Sec. 5 we present numerical estimations of the explicit solutions for the coupling. The cases of 
standard PT and of APT are separately considered. We compare numerically Pade and the iterative 
approximants for the three-loop coupling. The scope of validity for the iterative approximant is 
estimated. We give numerical results for the "analyticized couplings" A\(Q 2 , /) and 2li(s,/) to 
second and third orders. The differences between these quantities and a s (Q 2 , /), are estimated. The 
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numerical results for the global functions A n (Q 2 ) and 2l n (s) (n = 1, 2, 3), to second and third orders, 
(see Tables 7-12) are given. 



2 Exact solution for the two-loop coupling in the spacelike region 

The running coupling of QCD satisfies the differential equation [^ 

Q 2 &a fo2 f) = P f MQ 2 , /)) = -£ /«+ 2 (q 2 , /), (i) 

a s (fj, 2 ) = a s = |-, g is the renormalized coupling constant, is the renormalization point, and 

/ denotes the number of quark flavors. In the class of schemes where the beta-function is mass 

( f) ( f) f 

independent (3q and 0^ are universal and the result for f3 2 is available in the modified MS (MS) 

scheme 



/jZ-J-fll-M 0f-J_( 1O2 -™A 0f- -J- ( ?857 _ 5033/ 325£\ 
^° " 4vr V 3 V ' 1 ~ (4vr) 2 V 3 ^ ' ^ 2 ~ (4tt)3 ^ 2 18 + 54 J ' [ > 



For convenience in what follows we shall omit index f in the coefficients /3 k . Exact two-loop solution 



to Eq. <JT|) is given by 0, [12 1 



6 = /3i//3q, A = A-p^r and W(C) denotes the Lambert W function flj"7|] , the multivalued inverse of 
C = W"(C) exp W((). The branches of W are denoted Wk{C,),k = 0, ±1, ... . A detailed review of 
properties and applications of this special function can be found in [jT^ . The three-loop solution 
(with Pade transformed beta- function) for the coupling is [12| 



(3) fh2 f , ft 1 , 1 [Po^fQ 2 , 

Expressions (||) and (||) allows us to perform analytical continuation in the complex Q 2 plane and 
to calculate discontinuity along the negative Q 2 axis. [] In this way we construct the corresponding 
analytically improved expressions for the coupling. The "analyticized n-th power" of the coupling, 
obtained from the solution (||), can be written as 

aswj) = = ijf = i fc^^vnn*, ( 5) 

for < / < 6, the spectral function, p n (a, f) = prt, f) = Q{a s (— a — ^0)} n , is given by 

Z(t) = 1 exp^o^/ft 2 - t/*> + i(l/b - 1)tt). (7) 
oe 

Taking the limit /?2 — ^ in (|6|) one can reproduce the corresponding formula for the two loop case. 
b We use the notation Q 2 — —q 2 , Q 2 > corresponds to a spacelike momentum transfer. 



with 



c For details of analytical continuation we recommend papers [11-14]. 
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3 Continuation of the QCD running coupling to the timelike region 

The APT approach allows us to define the QCD coupling in the timelike region in a correct manner 



18, [L9(]. Here, instead of common power series for a timelike observables, there appears 



asymptotic series over the set of oscillating functions {2l n (s,/)} || 0. The functions 2l n (s, /) are 
defined by the elegant formula [ 20 1 



1 f°° da 

2ln(s,/) = -/ Pn((T,f), 

7T Js cr 



(8) 



here the spectral function is p n (a, f) = Q{a s (— a — z0)} n . In the one- loop order the set {2l n (s)} was 
studied in Ref.pl]. In this case, the first four functions of the set are given by 



I? 



PO 7T 



1 - 1 arctan(— )), 2lf } (s,/) ' 



1 



7T 



1 



In s 



/3 3 (ln 2 s + 7r 2 ) 2 ' 



4 [ ,J> /3 4 (ln 2 s + 7 r 2 ) 3 ' 



(9) 



here s = s/A 2 . In this section we will calculate 2l n (s, /) to second and third orders. Let us define 
the auxiliary functions 



Rn(sJ) 



1 

7T Jln i 



dta n (t,f). 



(10) 



where Qa n (t, f) = p n (t, f) = p n (cr, f) with a = exp(i). Then 2l n (s, /) = QR n (s, f). The expressions 
for 5, at the two and three loop orders, can be read from (0). In the two loop case 



: z(t) 



-l-t/b+itj> 



* = tt(J-1). 



fa 1 + Wi(z(t)) ■ b 
Integral (10), with (|TT|), can be rewritten as a contour integral in the complex z-plane 

rZs dz 1 



4 2 W) 



Pi 



here 



Pn 



■n-2 



on-1 ' 



1 



, z (l + WxCz))"' 



ee 



7T 



1 



1 



vr /3™ 

and the limit e — > is assumed. Let us introduce the new integration variable in (12) 

W(z) dz 



uj = W(z), 



doj 



then 



R n 2) {s,f)= Pn 
For n > 2, from (15) we find the relation 



l + W(z) z 
WxM l + UJ 



div. 



Pr n -^' J ' • (n-2) 
Eq. ( |l6| ) can be rewritten as the recurrence relation for 2l n (s) 



3(1 + Wi(z.)) 



2-n 



m ^ S s J) = " 2 )(Po%£-l( 8 > f) + WtiK*, /)), 



(11) 

(12) 
(13) 

(14) 
(15) 

(16) 
(17) 
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formula (|17| ) gives the generalization of the similar one-loop order relation obtained in paper 
From (|l7|) with the help of Eq. IM) we find analogical formula for the spectral function 



din a 



(2), 



(18) 



Let us multiply Eq. (18) by the factor (a + Q 2 ) 1 and take the integral over the region < a < oo. 
Integrating by parts and taking into account the condition p n _2(cr)cr/((3 2 + f)|JLo = we obtain 



MjSg(gVf) 

SlnQ 2 



-(n - 2)(/?o42i(Q 2 > /) + PiA^iQ 2 , /))• 



(19) 



Note that for n = 3 Eqs. (17) and ((T^) are analogous to the basic Eq. (|T]) with a™ replaced by 2l n and 
«4n respectively. We remark, that Eqs. (p~7|) - (|T9[) could also be derived on general grounds using the 
RG equation (|l]) together with the APT prescription. For this derivation there is no necessity in the 
explicit (exact) solutions of the RG equation. In higher orders, similar equations (which generalize 
Eq. ([[]) for A\ and 2li) can also be obtained [21]. In particular, to the k-th order, the equation holds 



dlnQ 2 



k-l 



-n 



(k) 

n+N+l 



(Q 2 J), 



(20) 



N=0 



the derivation of this equation will be given elsewhere [21 



(2) 

In the two-loop case, it is sufficient, to calculate 2li and 2I2 • Note that R\ (s,f) is divergent in 
the limit e — ► 0, (see ( |i2| ) ). Nevertheless, it has finite imaginary part By direct calculation we 
find 

((2)/ « 1 



npi \l + Wi{z s )J 



21 



(22) 



with z s given by 
s — > 0, we verify 



+ W 1 (z s ) / 

13|). Using the known asymptotic behavior of the W- function [17], in the limit 



a( 2 )( S ,/)=2lf } ( S ,/)^^, and ^\ S J)^0 for n > 1. 
Analogically in three loop case we find 

2li 3) (s, /) = a^(s, /) = — ^ (-%HWt(Z s )) + (1 - 1)9 ln(7/ + Wi(Z a )) - it) 

Wt(Z s ) 



4 3) (s,/) = ^r( 4$ in 



1 



7r/3i \n 



77 



1 

f? + ^(Z s 



(23) 



(24) 



(25) 



"Pc./)-^e.c./) +lK „_ a) 

where n > 2, 77 = 1 — P0P2/P1 and 



2-n 



^ -(i-^ + WiC^)) 1 -" (26) 



n — 1 ry 



Z s = I( s /A 2 )- 1 / fe exp(-^ + z(l/6- 1)tt). 



J for the asymptotic behaviour of the W function see paper 117 
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Note that the "analyticized" perturbative expansions for timelike observables (which contain specific 
functions 2l n ) may be rewritten as power series in traditional coupling a s (s) with modified by tt 2 - 
factors coefficients || [7J|. Previously, these modified power series have been obtained in [22, |23|| , 



Application of the series can be found in papers [23-29]. Thus, "7r 2 -effects" for various timelike 
quantities have been estimated in paper | |29y , in particular, it was found that the 7r 2 -factors give 
dominating contributions to the coefficients of R(s) = a to t(e + e~ — > hadrons)/o"(e + e~ — > fi + p,~). On 
the other hand, in recent paper |7j various timelike events was analyzed in the f=5 region. Higher- 
order "7r 2 -effects" have been taken into account properly. It was found that the extracted values for 
a s (M 2 ) are influenced significantly by these effects. 

4 Matching procedure and construction of the global space-like 
and time-like couplings 

In MS-like renormalization schemes important issue is how to introduce the matching conditions for 
the strong coupling constant at the heavy quarks thresholds. In literature few different recipes are 
known (see for example [||, |7], [H], [i^, [34|, |35|). Here, we follow works || Let us impose the 
continuity relations 

&s(Mf,Af-i, / — 1) = a s (Mj,Af, /). (27) 
Inserting in (27) the three- loop solution (||) we solve equation ( p7|) for Af 

A f = M f {-b f F(z f _x) exp(r ?/ + F(z f ^))} b f/ 2 ( 28 ) 
where b f = /?f/(/#) a , j = 1 - $ P f 2 / (p{) 2 , 

F{z f -i) = (Vf-l + W-iizf-ifi^tk ~ Vf- (30) 

Po Pi 

In paper || special model for the spectral functions Pn{&) was proposed 

Pn(a) = Pn (a, A 3 , 3) + ^e(a - M 2 )(p n (a, A f , f) - p n (a,A f ^, f- 1)), (31) 

here the mass Mt corresponds to the quark with flavor f, and Af is determined according formula 
(p8|). Inserting (|3l]) in formula (|8|) we find following expression for the "analyticized powers" of the 
global coupling in the timelike region Q 

2n(s) = an(«, f) + C n (f) for M f <yfs<M f+l , (32) 
where the shift coefficients C n {f) are defined by relation 

C n (f) = MMf+iJ + 1) - 2tn(M 2 +1 , /) + C k (f + 1) (33) 

with C n (6) = 0. The analogical formula follows for the corresponding global spacelike functions 

A n (Q 2 ) 

1 [ M l da 



. 1 f 4 da , A „, 

A n {Q 2 ) = - / ^ 2 p n {a, A 3 , 3 + 

1 A f M f+i da . , PN 1 Z* 00 da . A . 

_ 2^/ — Tia^A/,/ +- — T^Pn <t,A 6 ,6) 34 
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5 Numerical estimations 



For the quarks masses throughout this paper we assume the values Mi = M2 = M3 = 0, M4 = 1.3 
GeV, M 5 = 4.3 GeV and M 6 = 170 GeV. 

In practice, usually, the iterative approximation for the coupling [36] is used 



.(D,rf n_ 1 A lni , 1 , ft 



■w> = si - |if + ^ if 1 " 2 L - L - » + f J • (35) 

where L = mQ 2 / A fj s - The 

same formula is used in the timelike region. It is instructive to compare 

solutions (||) and ( |35| ) with exact numerical solution of the RG equation ctnm{Q 2 ,f)- Numerical 
results for these functions are summarized in the Table 2. We see that, for Q > 1 GeV, the expression 
(Q) is more accurate then ( |35|) and the difference between these functions becomes noticeable for 
Q < 1.2 GeV. 

In Table 3 we give results for the three- loop Pade approximated coupling a^ 3 \Q 2 ,5) and the 
corresponding analytic coupling A^\Q 2 , 5). In Table 4, a^(s,5) and St^ (s, 5) are compared. The 
interval 5 GeV < Q, 1/s < 200 GeV is chosen and it is assumed that A5 = 264 MeV. Here, we 
observe the inequality 

a^(Q 2 ,5) >A[ k) (Q 2 ,5) >%{ k) (Q 2 ,5), (36) 

the relative difference A (%) between a^ 3 \Q 2 ,5) and A^\Q 2 ,5) decreases from 1.4% at Q = 5 
GeV to 0.15% at Q = 200 GeV, whereas the difference between a^ 3 \s,5) and 21^ (s, 5) is more 
appreciable: A(%) = 7.5% at y/a = 5 GeV and A(%) = 1.6% at y/a = 200 GeV. The " mirror 
symmetry " || between ^4.^ (Q 2 ,5) and 21^ (s,5) is essentially violated (see Tables 3 and 4). Thus, 
the relative difference, A(%) = (^J 3) (Q 2 ,5) - 2tf \s, 5))/A^ \Q 2 , 5) * 100, monotonically decreases 
from 5.6 % at Q = yfe = 1 GeF to 1.5% at 200 GeF. 

In Table 5 various functions at the two- loop and three- loop orders are compared. We choose 
/ = 5 and A5 = 215 MeV . Then, from the matching formula (|28|), in the two-loop case (r) = 1) we 
find A3 = 363 MeV, while A3 = 340 MeV in the three-loop case. The interval 
20 GeV < Q, yfs < 170 GeV is chosen. Here, we demonstrate the stability of results of PT and of 
APT, with respect to the higher-loop corrections. In this region, the numerical difference between 
the three-loop and the corresponding two- loop couplings are of order 0.3% — 0.2%. 

In Tables 7-12 we have summarized numerical results for the global three- loop functions , 
An\Q 2 ) and %in\s) for n=l,2,3. The values for the parameter A3 are chosen at 350 MeV, 400 
MeV, and 450 MeV. The corresponding values for the Af for / > 3 are calculated from the 
three-loop matching condition (^), (see Table 1). 

Let us compare the global three-loop function A^{\Q 2 ) (see Table 9) and a^ 3 \Q 2 ,5) (Table 3) 
in the case A3 = 400 MeV (the corresponding value for A5 is 264 MeV). We see that the global 
coupling A^{\Q 2 ) does not obey the inequality (j36|). The difference a^(Q 2 ,5) — A^(Q 2 ) becomes 
negative when Q increases. For A3 = 400 MeV, this takes place at Q ~ 17 GeV. This difference 
is small but it increases with Q. It is about 0.6% at Q = 200 GeV. The same effect is occurred for 
other values of the A3. This enhancement of the global APT coupling can be explained from formula 
(|34]). It is easy to verify that (|34| ) contains additional positive non-perturbative contribution, which 
was not occurred in the case of the local APT coupling (|5|). This contribution increase the coupling 
for large values of Q. 

(3) 

Such a behavior does not occur in the case of the global timelike coupling W n (s) (compare 
Tables 4 and 10). However, the difference between a^(s,f = 5) and 21^ (s) is more appreciable, it 
is about 7% at s = 5 GeV, 3.4% at s = 20 GeV and 1.75% at s = 90 GeV. It is about 1% at s = 200 
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GeV. The relative difference between Af ] (Q 2 ) and 2tf \s) decreases from 6.9% at Q = yfs = 2 GeV 
to 2% at 100 GeV (see Tables 9 and 10). The relative differences A n (%) between An\Q 2 ) and 
2l^ 3) (s), for n=2 and 3, are even large. Thus, A 2 (%)=6.1% and A 3 (%)=11.1% at Q = y/s = 100 
GeF. 

With the algebraic computer system Maple V release 5 we were able to calculate the quantities 
{a s (Q 2 , f)} n and 2l n (s,/) (see formulas ([|),(||),(|4|)-(|26|)) with any arbitrary given accuracy. How- 
ever, this is no case for formulas (||) and (0). These integrals are singular at t — > dboo, therefore, 
one needs to use a cutoff. With Maple V release 5 the cutoff may be as big as 10 4 . This guarantees 
to obtain 4-5 reliable digits after decimal point. Most of our calculations we have performed with 
this precision. To obtain more accurate results we suggest following formula 

A n (Q 2 , f) = A n (Q 2 , f, R) + 2l n (A 2 e R , /) + 0(e~ R ), (37) 

Here A n (Q 2 ,f,R) denotes the regulated integral @: the integral is taken over the finite region 



R < t < R. We remark that formula (37) provides sufficiently high precision even for low values 



of the cutoff. In addition, for practical calculations, here we suggest the formula 

/o (s + Q 



A n (Q J) = Q I - I ^ 2 , 2 gn(s,/), (38) 



evidently, this relation is also valid for the global quantities, A n (Q ) and 2t„(s 



6 Conclusion 

The "analyticized powers" of the coupling in the spacelike region A^\Q 2 ,f), to second and third 
orders, are written in terms of the Lambert W function (see integrals (||) and (|3§|)). 

The "analyticized powers" of the coupling in the timelike region, 21^ (s,f) (k=2,3), are ana- 
lytically calculated in terms of the Lambert W function (see formulas (|2l])-(|2^)). The recurrence 
relations for 21^ (s,/) and An\Q 2 , f) are derived (see formulas ( |i~7| ) and (Jl9[) ) . 

The matching conditions for crossing the quark flavor thresholds (|27| ) are solved explicitly for A f 
in the cases of the exact 2- loop and Pade improved 3-loop solutions (see formulas @, (^) and (|28|)). 

The global model for the coupling, of Refs.p-0, is considered up to the third order in the 
context of obtained explicit solutions. 

In Sec. 5 numerical estimations of the explicit solutions for the powers of the standard and ana- 
lytical couplings are given. We have compared various solutions in the large region of momentum 
transfer and energy, 1 GeV < yfs,Q < 200 GeV (see Tables 1-12). We have confirmed that the 
differences between the powers of the standard iterative solution (S3) and the "analyticized timelike 
powers" 2ln (s) are appreciable even for moderate energies. 
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Table 1: The values of Af determined by the threshold matching condition (|28|). A3 is chosen as a 
basis quantity. 



The three loop case 


Exact two loop case 


A3 


A 4 


A 5 


A 6 


A3 


A 4 


A 5 


A 6 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


300 


258.4 


183.0 


76.2 


300 


248.4 


170.4 


69.2 


310 


268.2 


190.8 


79.7 


310 


257.6 


177.3 


72.3 


320 


278.0 


198.6 


83.4 


320 


266.8 


184.3 


75.5 


330 


287.9 


206.5 


87.0 


330 


276.0 


191.4 


78.7 


340 


297.8 


214.4 


90.8 


340 


285.3 


198.5 


81.9 


350 


307.8 


222.5 


94.6 


350 


294.6 


205.7 


85.2 


360 


318.0 


230.6 


98.4 


360 


303.9 


212.9 


88.5 


370 


328.0 


238.8 


102.3 


370 


313.2 


220.1 


91.9 


380 


338.1 


247.1 


106.2 


380 


322.6 


227.4 


95.2 


390 


348.4 


255.5 


110.2 


390 


332.0 


234.8 


98.7 


400 


358.7 


263.9 


114.3 


400 


341.4 


242.2 


102.1 


410 


369.0 


272.5 


118.4 


410 


350.9 


249.6 


105.6 


420 


379.4 


281.1 


122.6 


420 


360.4 


257.1 


109.1 


430 


389.9 


289.8 


126.8 


430 


369.9 


264.7 


112.7 


440 


400.5 


298.6 


131.1 


440 


379.5 


272.2 


116.2 


450 


411.1 


307.5 


135.4 


450 


389.0 


279.9 


119.8 


460 


421.8 


316.5 


139.8 


460 


398.6 


287.5 


123.5 


470 


432.6 


325.6 


144.3 


470 


408.3 


295.2 


127.2 


480 


443.4 


334.8 


148.8 


480 


417.9 


303.0 


130.9 


490 


454.4 


344.1 


153.4 


490 


427.6 


310.8 


134.6 


500 


465.4 


353.5 


158.1 


500 


437.3 


318.6 


138.4 
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Table 2: The three loop Pade, a d !, and the iterative, aq t , couplings versus the exact numerical 

(3) 

solution of the RG equation, c4m- We take A3 = 400 MeV and the values of Af (/ > 3) are given 
in Table 1. We denote A pd (%) = \o£l - c$\/ofil * 100, , A it (%) = \o£l - c$\/c£l * 100 and 

a p ,(%) = i«S ) -4 3) i/«S*ioo 



Q Gel/ 


„ (3) m2\ 


„ (3) /,n2\ 


„ (3) /,o2\ 


A 




A 


GeV 


f=3 


f=3 


f=3 


{ fry \ 

(%) 


(%) 


(%) 


.8 


.76491 


.90931 


.88340 


18.9 


15.5 


2.9 


.9 


.63323 


.68022 


.69179 


7.4 


9.3 


1.7 


1.0 


.55414 


.57854 


.58784 


4.4 


6.1 


1.6 


1.1 


.50028 


.51548 


.52195 


3.0 


4.3 


1.3 


1.2 


.46075 


.47126 


.47589 


2.3 


3.3 


1. 


1.3 


.43025 


.43803 


.44153 


1.8 


2.6 


.8 


1.4 


.40587 


.41191 


.41469 


1.5 


2.2 


.7 


1.5 


ooroo 

.38583 


.39069 


.39301 


1.3 


1.9 


.6 


1.6 


.36901 


.37302 


.37503 


1.1 


-1.6 


.5 


1.7 


.35464 


.35803 


.35982 


1. 


1.5 


.5 


1 8 


.34220 


.34512 


34674 


9 


1 3 

J. . O 


. O 


1.9 


.33130 


.33384 


33535 


.8 


1.2 


.5 


2 


32165 

. O — 1 \) O 


32389 

. O t.) (J .J 


32530 

. '.) — • O O V 7 


7 


1.1 


.4 


2.1 


.31303 


.31503 


.31637 


.6 


1.1 


.4 


2.2 


.30527 


.30707 


30835 

. OOOOO 


.6 


1.0 


.4 


2.3 


.29824 


29988 


.30111 


.6 


1. 


.4 


2.4 


.29184 


.29334 


.29452 


.5 


.9 


.4 


2.5 


28598 


.28735 


28850 


.5 


.9 


.4 


2.6 


.28058 


.28185 


.28297 


.5 


.9 


.4 




a ( nl(Q 2 ) 


a {3 ](Q 2 ) 


cv ( f } (Q 2 ) 




A it 




GeV 


f=4 


f=4 


f=4 






(%) 


2 


33260 


33392 


33392 

• OOO l7 4J 


.4 


.4 





3 


.27479 


.27540 


.27565 


.2 


.3 


.1 


4 


.24527 


.24565 


.24598 


.2 


.3 


.1 


5 


.22664 


.22691 


.22726 


.1 


.3 


.2 


6 


.21350 


.21372 


.21406 


.1 


.3 


.2 


7 


.20359 


.20376 


.20410 


.1 


.3 


.2 


8 


.19575 


.19590 


.19623 


.1 


.3 


.2 


9 


.18935 


.18948 


.18980 


.1 


.2 


.2 


10 


.18398 


.18410 


.18441 


.1 


.2 


.2 


Q GeV 


a£i(Q 2 ) 




*fiQ 2 ) 




A lt 




GeV 


f=5 


f=5 


f=5 


(%) 


(%) 


(%) 


10 


.18845 


.18849 


.18847 


.02 


.01 


.01 


15 


.17126 


.17128 


.17130 


.01 


.03 


.01 


20 


.16090 


.16092 


.16096 


.01 


.03 


.02 


25 


.15372 


.15373 


.15378 


.01 


.04 


.03 


30 


.14832 


.14833 


.14838 


.01 


.04 


.03 
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Table 3: The Q 2 dependence of the three loop Pade improved coupling a^ 3 \Q 2 ,f), (f§), and 
corresponding analytic coupling Af\Q 2 , f) for / = 5 and A 5 = 264 MeV (A3 = 400 MeV). A 
denotes the relative difference between the couplings. 
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Table 4: The three loop analytic coupling 21^ (s, 5) ( see (pi])) versus the ordinary 3-loop Pade 
approximated coupling a^(s, 5). We have assumed that A3 = 400 MeV, correspondingly A5 = 264 
MeV. 
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Table 5: The two and three loop couplings a^(Q 2 ,5), A { ^ } (Q 2 ,5) and 2lj fe) (,s,5) (A; = 2,3) as a 
functions of variables Q and \fs. We take f=5 and A 5 = 215 MeV. 
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.11483 


.11298 


.11523 


.11506 


.11319 


115 


.11425 


.11408 


.11228 


.11448 


.11431 


.11248 


120 


.11355 


.11338 


.11161 


.11377 


.11360 


.11181 


125 


.11288 


.11271 


.11098 


.11310 


.11293 


.11117 


130 


.11225 


.11208 


.11037 


.11246 


.11230 


.11057 


135 


.11164 


.11148 


.10980 


.11186 


.11169 


.10999 


140 


.11107 


.11090 


.10926 


.11128 


.11111 


.10945 


145 


.11052 


.11035 


.10873 


.11073 


.11056 


.10892 


150 


.10999 


.10983 


.10823 


.11020 


.11003 


.10842 


155 


.10949 


.10932 


.10776 


.10969 


.10953 


.10794 


160 


.10901 


.10884 


.10730 


.10921 


.10904 


.10748 


165 


.10854 


.10838 


.10685 


.10874 


.10857 


.10703 


170 


.10810 


.10793 


.10643 


.10829 


.10813 


.10660 
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Table 6: The shift constants Ck{f) as a functions of the parameter A3. 



A 3 MeV 


250 


300 


350 


400 


450 


d(3) 


.0110 


.0137 


.0169 


.0203 


.024 


C 2 (3) 


.0062 


.0079 


.0099 


.0120 


.014 


Ca(3) 


.0022 


.0028 


.0035 


.0042 


.0049 


d(4) 


.0026 


.0032 


.0037 


.0043 


.0049 


C 2 (4) 


.0013 


.0016 


.0019 


.0023 


.0027 


C 3 (4) 


.0004 


.0005 


.0007 


.0008 


.0010 


d(5) 


.0003 


.0003 


.0003 


.0083 


.0004 


C 2 (5) 


.0001 


.0001 


.0001 


.0020 


.0001 


Cs(5) 


.0000 


.0000 


.0000 


.0003 


.0000 
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Table 7: The three-loop Pade approximated global analytic coupling A\{Q 2 ) and the corresponding 
global "analyticized powers" A2{Q 2 ) and A%{Q 2 ) as a functions of the momentum transfer Q for 
A 3 = 350 MeV . The matching conditions give A 4 = 307.8 MeV, A 5 = 222.5 MeV and A 6 = 94.6 
MeV. 



Q GeV 


Ai{Q 2 ) 


A 2 (Q 2 ) 


A 3 (Q 2 ) 


Q GeV 


A^Q 2 ) 


A 2 (Q 2 ) 


A 3 (Q 2 ) 


1 


.36571 


.09320 


.015540 


38 


.13774 


.01898 


.00258 


2 


.28689 


.06861 


.012776 


40 


.13650 


.01864 


.00252 


3 


.25144 


.05641 


.010598 


42 


.13534 


.01833 


.00246 


4 


.23041 


.04903 


.009107 


44 


.13442 


.018047 


.00240 


5 


.21615 


.04403 


.008042 


46 


.13340 


.01777 


.00235 


6 


.20566 


.04039 


.007245 


48 


.13244 


.01752 


.00230 


7 


.19753 


.037601 


.006626 


50 


.13153 


.01728 


.00225 


8 


.19098 


.035378 


.006131 


52 


.13067 


.01706 


.00221 


9 


.18555 


.033557 


.005724 


54 


.12986 


.01685 


.00217 


10 


.18095 


.032033 


.005384 


56 


.12908 


.01665 


.00213 


11 


.17698 


.030734 


.005095 


58 


.12834 


.01646 


.00210 


12 


.17351 


.029610 


.004845 


60 


.12763 


.01628 


.00206 


13 


.17043 


.028626 


.004628 


62 


.12696 


.01611 


.00203 


14 


.16768 


.027755 


.004436 


64 


.12631 


.01595 


.00200 


15 


.16520 


.026977 


.004266 


66 


.12569 


.01580 


.00197 


16 


.16295 


.026277 


.004113 


68 


.12510 


.01565 


.00195 


17 


.16089 


.025643 


.003975 


70 


.12452 


.01551 


.00192 


18 


.15900 


.050652 


.003851 


72 


.12397 


.01537 


.00190 


19 


.15725 


.024535 


.003737 


74 


.12344 


.01524 


.00187 


20 


.15563 


.02405 


.00363 


76 


.12293 


.01511 


.00185 


21 


.15412 


.02359 


.00353 


78 


.12244 


.01499 


.00183 


22 


.15270 


.02317 


.00344 


80 


.12196 


.01488 


.00181 


23 


.15138 


.02278 


.00336 


82 


.12150 


.01477 


.00179 


24 


.15013 


.02242 


.00328 


84 


.12105 


.01466 


.00177 


25 


.14895 


.02208 


.00321 


86 


.12062 


.01455 


.00175 


26 


.14784 


.02176 


.00315 


88 


.12020 


.01445 


.00173 


27 


.14678 


.02145 


.00308 


90 


.11979 


.01436 


.00171 


28 


.14578 


.02117 


.00303 


92 


.11940 


.01426 


.00170 


29 


.14483 


.02090 


.00297 


94 


.11902 


.01417 


.00168 


30 


.14392 


.02064 


.00292 


96 


.11864 


.01408 


.00166 


32 


.14222 


.02016 


.00282 


98 


.11828 


.01400 


.00165 


34 


.14066 


.01973 


.00273 


100 


.11793 


.01391 


.00163 


36 


.13923 


.01934 


.00265 
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Table 8: The three-loop Pade approximated global analytic coupling 2li(s) and the corresponding 
global "analyticized powers" ^(s) and ^(s) as a functions of the energy y/s for A 3 = 350 MeV. 



fa Cf>V 


91, f <A 

-ai^sj 


a} 
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a\ 
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40 


.13326 
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.00220 


165 


.10785 


.01148 


.00120 


42 


.13219 


.01708 


.00215 


170 


.10742 


.01139 


.00119 


44 


.13118 


.01682 


.00210 


175 


.10704 


.01131 


.00118 


46 


.13023 


.01659 


.00206 


180 


.10667 


.01124 


.00116 


48 


.12933 


.01637 


.00202 


185 


.10632 


.01116 


.00115 


50 


.12849 


.01616 


.00198 


190 


.10598 


.01109 


.00114 


52 


.12768 


.01596 


.00195 


195 


.10565 


.01102 


.00113 


54 


.12692 


.01578 


.00192 


200 


.10533 


.01096 


.00112 
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Table 9: The three-loop Pade approximated global analytic coupling Ai(Q 2 ) and the corresponding 
"analyticized powers" ^(Q 2 ) and As(Q 2 ) as a functions of the momentum transfer Q for A3 = 400 
MeV. The matching conditions give A 4 = 358.7 MeV and A 5 = 263.9 MeV. 



Q GeV 


A^Q 2 ) 


A 2 (Q 2 ) 


MQ 2 ) 


Q GeV 


AxiQ 2 ) 


A 2 (Q 2 ) 


A 3 (Q 2 ) 


1 


.38656 


.09988 


.01629 


58 


.13190 


.01743 


.00228 


2 


.30263 


.07438 


.01387 


60 


.13116 


.01723 


.00224 


3 


.26444 


.06128 


.01164 


62 


.13044 


.01705 


.00221 


4 


.24171 


.05325 


.01005 


64 


.12976 


.01687 


.00217 


5 


.22627 


.04777 


.00890 


66 


.12910 


.01670 


.00214 


6 


.21492 


.04377 


.00803 


68 


.12847 


.01654 


.00211 


7 


.20613 


.04069 


.00734 


70 


.12786 


.01639 


.00208 


8 


.19905 


.03824 


.00679 


72 


.12728 


.01624 


.00206 


9 


.19318 


.03623 


.00634 


74 


.12672 


.01610 


.00203 


10 


.18822 


.03454 


.00596 


76 


.12618 


.01596 


.00200 


11 


.18394 


.03311 


.00564 


78 


.12566 


.01583 


.00198 


12 


.18020 


.03187 


.00536 


80 


.12515 


.01570 


.00196 


13 


.17689 


.03078 


.00490 


82 


.12467 


.01558 


.00193 


14 


.17393 


.02982 


.00471 


84 


.12420 


.01547 


.00191 


15 


.17127 


.02896 


.00454 


86 


.12374 


.01535 


.00189 


16 


.16885 


.02819 


.00438 


88 


.12330 


.01524 


.00187 


17 


.16664 


.02749 


.00424 


90 


.12287 


.01514 


.00185 


18 


.16461 


.02685 


.00411 


92 


.12245 


.01504 


.00183 


19 


.16273 


.02627 


.00400 


94 


.12204 


.01494 


.00182 


20 


.16099 


.02573 


.00389 


96 


.12165 


.01484 


.00180 


21 


.15937 


.02523 


.00379 


98 


.12127 


.01475 


.00178 


22 


.15786 


.02477 


.00370 


100 


.12090 


.01466 


.00177 


23 


.15644 


.02434 


.00361 


105 


.12001 


.01445 


.00173 


24 


.15510 


.02394 


.00353 


110 


.11917 


.01425 


.00169 


25 


.15385 


.02357 


.03460 


115 


.11839 


.01406 


.00166 


26 


.15266 


.02321 


.00339 


120 


.11765 


.01389 


.00163 


27 


.15153 


.02288 


.00332 


125 


.11695 


.01372 


.00160 


28 


.15046 


.02257 


.00326 


130 


.11628 


.01357 


.00157 


29 


.14944 


.02227 


.00320 


135 


.11565 


.01342 


.00155 


30 


.14847 


.02199 


.00309 


140 


.11505 


.01328 


.00153 


32 


.14666 


.02147 


.20099 


145 


.11448 


.01315 


.00150 


34 


.14500 


.02100 


.00290 


150 


.11394 


.01303 


.00148 


36 


.14347 


.02056 


.00282 


155 


.11341 


.01291 


.00146 


38 


.14205 


.02017 


.00275 


160 


.11291 


.01280 


.00144 


40 


.14074 


.01980 


.00269 


165 


.11243 


.01269 


.00142 


42 


.13951 


.01946 


.00262 


170 


.11197 


.01258 


.00141 


44 


.13835 


.01915 


.00256 


175 


.11153 


.01248 


.00139 


46 


.13727 


.01885 


.00251 


180 


.11110 


.01239 


.00137 


48 


.13625 


.01858 


.00245 
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Table 10: The three-loop Pade approximated global analytic coupling 2li(s) and the corresponding 
global "analyticized powers" ^(s) and ^(s) as a functions of the energy y/s for A 3 = 400 MeV. 
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Table 11: The three-loop Pade approximated global analytic coupling A\{Q 2 ) and the corresponding 
global "analyticized powers" A2{Q 2 ) and A%{Q 2 ) as a functions of the momentum transfer Q for 
A 3 = 450 MeV. The matching conditions give: A 4 = 411.1 MeV, A 5 = 307.5 MeV and A 6 = 135.4 
MeV. 
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Table 12: The three-loop Pade approximated global analytic coupling 2li(s) and the correspond 
"analyticized powers" 2l2(s) and 2ts(s) as a functions of the energy y/s for A 3 = 450 MeV . 
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